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Making use of the complete calculation [1] of the chiral six-point correlation function 

C{z) = {(/il,2(?!)l,2"I>l/2,o(2,2)<^l,2<^l,2) , 



with the four 01,2 operators at the corners of an arbitrary rectangle and the point z = a; + iy in 
0^ , the interior, for arbitrary central charge (equivalently, SLE parameter k > 0), we calculate various 

04 ■ quantities of interest for percolation (k = 6) and many other two-dimensional critical points. In 

particular, we use C to specify the density at z of critical clusters conditioned to touch either or 
' (— I ' , both vertical sides of the rectangle, with these sides 'wired,' i.e. constrained to be in a single cluster, 

^ • and the horizontal sides free. These quantities probe the structure of various cluster configurations, 

I ' including those that contribute to the crossing probability. 

, We first examine the effects of boundary conditions on C for the critical 0(n) loop models in both 

^ ■ high and low density phases and for both Fortuin-Kasteleyn (FK) and spin clusters in the critical 

Q-state Potts models. A Coulomb gas analysis then allows us to calculate the cluster densities 
with various conditionings in terms of the conformal blocks calculated in [1]. Explicit formulas 
generalizing Cardy's horizontal crossing probability to these models (using previously known results) 
are also presented. 

These solutions are employed to generalize previous results demonstrating factorization of higher- 
order correlation functions to the critical systems mentioned. An explicit formula for the density 
' of critical percolation clusters that cross a rectangle horizontally with free boundary conditions is 

\^ ' also given. Simplifications of the hypergeometric functions in our solutions for various models are 

l/^ , presented. 

High precision simulations verify these predictions for percolation and for the Q = 2 and 3-state 
Potts models, including both FK and spin clusters. Our formula for the density of crossing clusters 
in percolation in open systems is also verified. 
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I. INTRODUCTION 

The conformal invariance of crossing probabilities in critical percolation was first suggested by Aizenman, on the 
basis of numerical results by Langlands et al. [2]. This prompted Cardy [3] to apply the methods of boundary 
conformal field theory (CFT) to derive his celebrated formula for the horizontal crossing probability in a rectangle of 
arbitrary aspect ratio. This formula is an excellent example of the predictive power of CFT. 

Recently there has been a renewed interest in crossing probabilities as prototype non-local observables; rigorous 
proofs have appeared [4] , and progress has been made in exploiting CFT and other methods to capture non-local phe- 
nomena [5, 6]. This includes our previous crossing results for percolation [7]. Much of the contemporary research has 
been motivated by the development of Schramm-Loewner Evolution (SLE) [8], which provides a rigorous framework 
for conformal symmetry but with a perspective contrasting CFT. Work on crossing probabilities has also stimulated 
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new research in number theory, in particular the development of the theory of higher-order modular forms. See [9] 
and [10] for details and references. 

In [11] and [12] we established, by use of conformal field theory and high-precision simulation, exact and universal 
factorizations of certain higher-order correlation functions in terms of lower-order correlation functions for percolation 
clusters in two dimensions at the percolation point. In that work, the correlation functions involved the density of 
critical percolation clusters constrained to touch one or two isolated boundary points, or single boundary intervals. 
The question of which conformal correlation functions factorize in this way is examined in [13]. 

A more recent paper [14] considers percolation densities in a rectangle conditioned to touch one or both vertical 
sides, with the sides 'wired', i.e. constrained to belong to a single cluster. Evaluating these quantities via CFT requires 
calculation of the correlation function (1) 

C{Z) = (</'t2(O)0?,2(i)'i>l/2,o(^,5)'/'t2(^)'/'t2(^ + i)>7^ , 

in a rectangular geometry TZ with the boundary operators at the corners. This is tantamount to the challenging 
task of evaluating a chiral six-point function in the complex plane. In [1] all solutions of the PDEs governing this 
correlation function are determined, for arbitrary central charge (or equivalently, SLE parameter k > 0) in terms of 
algebraic and Appell hypergeometric functions. 

The solutions found in [1] are employed here to give weights and cluster densities for various quantities in a range 
of critical models. (We use the term "cluster density" in a particular way; the definition is given in subsection II B.) 
We index the models via the SLE parameter k in this article, as it includes both branches of the 0(n) model and thus 
the FK and Potts spin cluster models. The interpretation of $1/2,0 ^ density operator holds as long as its weight 
is positive, or when 8/3 < k < 8 (corresponding to < n < 2). Outside of this region, we do not offer a physical 
interpretation for C(z), though the solutions found in [1] remain valid. 

We also employ the results for C to obtain generalized factorization of correlations. In [14], which treats percolation 
only, a certain universal ratio of cluster densities within a rectangular region is considered, involving the densities of 
clusters touching the left and/or right sides of the rectangle with wired, or fixed, boundary conditions on those sides. 
The ratio is nearly constant everywhere in the rectangle, varying by less than 3%, and in addition does not depend 
on the vertical position y in the rectangle. In this paper, we generalize the percolation results to the range of critical 
models mentioned. How exactly the factorization holds depends on the model, but the independence from y is always 
valid. 

From a broader perspective, C(z) provides, via the densities, information on the structure of crossing and related 
cluster configurations (or equivalently loop configurations) in the critical models mentioned. In this sense it generalizes 
and deepens Cardy's analysis of crossing in percolation [3]. 

In section II we recall the results for the correlation function (1) obtained in [1]. The interesting independence from 
one coordinate mentioned (see (21)) appears here. Subsection II B then examines questions that arise in applying our 
results to the specific critical models mentioned. In particular, we examine the effects of boundary conditions, and 
obtain formulas for properly normalized densities. In addition, we give explicit formulas for a certain generalization 
of Cardy's horizontal percolation crossing probability [3]. Subsection II C identifies the solutions for the densities of 
various cluster configurations of interest by use of Coulomb gas methods. 

In section III we analyze a certain universal ratio p of densities, which (in the scaling limit) is equal to a ratio 
of involving critical cluster densities and crossing probabilities. This quantity is universal, and almost constant 
everywhere in the rectangle, regardless of aspect ratio, in many models, which implies an almost exact factorization of 
the corresponding correlation functions (or densities). This work generalizes the ratio examined in [14] for percolation. 

Section IV returns to the case of percolation. Here we take advantage of the "locality" property of percolation, as 
it is referred to in SLE (equivalent to the independence of local sites or bonds in lattice models) to derive a formula 
for the density of horizontal crossing clusters Pa^ in rectangles with free boundary conditions on all sides. 

In section V we compare our predictions for the ratio p, which governs factorization, and our formula for the density 
of percolation clusters crossing a rectangle with high-precision simulations. The agreement is very good. 

Section VI contains a detailed summary of our results. 

The correlation function C'(z) determines various densities; a similar correlation function without the <&i/2,o deter- 
mines crossing weights. These in turn depend on the Appell or hypergeometric functions, respectively. Appendix 
A presents results for various critical models where the hypergeometric functions simplify. 
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II. THEORY 
A. Results for the correlation function 

In this subsection we recall results obtained in [1] for the six-point correlator 

C{z) = {ct^lM4>l2(^)^i,2.,o{z, -zWi,2{RWi,2{R + i))7^ , (1) 
in the rectangular geometry TZ :~ {z ~ x + ly ^ <C\Q < x < R,Q < y < 1}. The aspect ratio is given by 

_ K{m) 

where K'{m) := K{l — m), with K is the complete elliptic integral. Conversely, the elliptic parameter m specifies the 
aspect ratio R via 

(Note that m differs from the standard modular lambda parameter, which is 1 — m here.) 
The conformal dimensions and central charge are 

(8-^)(3^-8) 

ftl/2,0 - /il/2,0 - ^ (4j 

h,.2 = ^ (5) 

K,, = ^ (6) 

K 

^ (3/.-8)(6-k) 

2k ' ^ ' 

where k is the Schramm-Locwncr Evolution (SLE) parameter and hi ^, is used below. Making use of the coordinates 

^ -.^ iin{xK'\mf , and := sn{y K'\l — mf , (8) 

we find that any solution that is a single conformal block can be written either in the form 

C(z) = /(e,^,m)G(C,m) , (9) 

where the algebraic prefactor / is given by 



f{^,-ip,m) = c(m) 
with c(m) given by 



(1-771^2)2 (1- (l-m)V'2)2 ^ 1 -''i.3/2+'^i/2,0 



$(l-0(l-mO 7/'(l-V^)(l-(l-m)V) 



(10) 



c(m) = 2'"^'^ (x')8''i.2+2''i/2,o (^(1 _ m)fh,,^ ^ (^^^ 



or as in (9) with G {£,, m) — > G (i/^, 1 — to). 

The coefficient 2''^-^ indicates the presence of a corner operator. The prefactor / is a K-depcndent power of the 
density of clusters in a rectangle attached to a fixed boundary. See [1] for details. 

The algebraic prefactor / is independent of boundary conditions. The conformal block G, on the other hand, is 
strongly dependent on them, as discussed below. This plays an important role in our results. Now G is given by one 
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of the five solutions 

r(2 - 8/k)T{16/k. - 1) [m(l - m)f^'' ^8/«-i/2 



Gi(e,m 
Gii(^, m 
Gm{^, m 



r(i2/.c)r(i - 4/«) [(1 _ _ 



Fi 1 



4 4 4 12 



(1 -m)^/'^ 

m6/«-i[e(l- 0(1-^6] 
(1 - mOi^/""3/2 

r(2 - 8/K)r(i6/K- 1) 



4/^^^i ;2-- 



1 — 771, 1 — m^^ , 
m, 



Fi 1 



4 4 



^ 16 „ 
,2-— ;2- 



1 — TO, 



2/''(X _ 08/K-1/2 



r(12/K)r(l - 4/k) (1 - m)6/«-i^4/K-i/2(i _ 
where Fi is the first Appell function 



1 — TO 

1 - toC 
F (1 4.4 4.12 

^'^ \ K K K K 



and 

m(l-0 1-C 



1 — TO^ ' 1 — TO^ 



Fi{a;bi,b2;c\zi,Z2) := ^ 

ij=0 



{a)i+j{bi)i{b2)j zi^ Z2^ 



i+j 



(12) 
(13) 
(14) 
(15) 

)■ 

(16) 
(17) 



with the Pochhammer symbol (z)„ := T{z + ri)/r(z). 

For our ranges of ^ and to values (0 < ^to < 1), (12-16) exhaust the convergent Frobenius series solutions to the 
differential equations that can be expressed with a single Fi. One can also use one of these solutions (or the five 
mentioned below) with — ^ t/;, to — > 1 — m}. Each of these is a single conformal block. Only three of them are 
independent, as discussed in section II A of [1]. 

There are also five other convergent solutions that can be expressed with a single second Appell function F2 (see 
[1] for details), of which we use only one here: 



Gyi{£_, to) := Gii - nGi = Giv — nGy 

r(2 - 8/K)r(i6/K - i)r(4/K) 



(18) 



to(1 - m)]^/" 



^(1 _ ^)-|8/'=-i/2 



1 — TO^ 



Fo 



r(i -4//^)r(8/K)2 

where n is the parameter of the 0{n) loop models, given by 

n = — 2 cos(47r/K) , 

with F2 defined by 

{a)i+j{bi)i{b2)jZi''z2^ 



16 ^.448 



^(1 - to) 

1-me 



(19) 



F2 (a;&i,62;ci,C2|zi,Z2) := ^ 



i!j!(ci),(c2)j 



(20) 



The fact that these solutions for G only depend on two variables rather than three follows on combining the 
conformal PDEs. One finds that 



(21) 



This equation is the basis for the y-independence of the factorization behavior discussed below. It indicates the 
presence of an unknown symmetry. 

For use below, we also recall the behavior of ^, ^ and to under the symmetries of the rectangle: mirror symmetries 
about X = R/2, y = 1/2 and x = y. These symmetry operations translate, respectively, into 



1-e 
1 - TO^ 



1 - V 

1 — (1 — to)?/) ' 



m] ,{1^,^,1-171) 



(22) 
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B. Physical interpretation of equation (1) and boundary conditions 

The correlation function (1) describes cluster densities for a variety of critical statistical mechanics models; in this 
section we examine this correspondence. We focus in particular on the Potts models and their boundary conditions 
in order to obtain explicit formulas. These results allow us to implement simulations (see section V) to test our 
predictions. We begin with critical 0(n) loop models (in the continuum limit) because they furnish a clear picture 
for these densities, and describe a continuum of critical points with < n < 2. There are many realizations of these 
models but in each case the degrees of freedom are closed fractal loops, each contributing to the partition function 
with a fugacity n. There are two critical branches depending on the energy cost per loop length, representing dilute 
and dense loop phases. As noted, we use the SLE parameter k to index our models. This is convenient as k nicely 
parameterizes both critical branches of the 0(n) loop model; the two parameters are related by (19) with k < 4 (resp. 
K > 4) corresponding to the dilute (resp. dense) phase. 

We can condition the loop ensembles so that open loop segments emerge from pairs of points on the boundary. 
We call these loops boundary arcs to distinguish them from the closed loops in the bulk. These boundary arcs are 
equivalent to SLE traces, and in CFT are implemented by placing 0i^2 Kac operators at their endpoints on the 
boundary. Because they are distinguished from bulk loops, we can give the boundary arcs a fixed weight of 1. The 
correlation functions we consider correspond to an 0(n) model with boundary arcs attached to the four corners of an 
arbitrary rectangle. 

In the 0(n) model, bulk points are said to be adjacent to a given boundary segment whenever a path can be drawn 
from the point to the segment without crossing any loops or boundary arcs. It has been argued using Coulomb gas 
techniques that the density operator for these adjacent points is $1/2,0(^1 ^) [15]. Mathematically our conformal block 
expressions ((12)-(16) and (18)) are valid for all k > 0, but we only discuss physical interpretations for 8/3 < k < 8. 
Outside of this range the density operator in the 0(n) models has weight zero, while <I'i/2,o does not. 

In order to proceed, we also need to consider generalizations of Cardy's crossing probability. They provide an 
example of applying boundary conditions appropriate to the Potts models, and a necessary normalization. We discuss 
these matters in subsections II B 1 and II B 2 below. 

In applying our results to model systems, three kinds of quantities arise. Unnormalized weights, which are denoted 
by n, crossing probabilities denoted by P, and cluster densities, also denoted by P. The crossing probabilities are true 
probabilities, and can be completely determined from our CFT results. For a given model, they depend on the aspect 
ratio R only. A cluster density (or more simply, "density" ) is a bit more complicated. By "density" we understand 
the scaled probability in a small region around a point with some specified conditioning, e.g. connected to a specific 
boundary segment. This is the probability, on a lattice with mesh size S, that a neighborhood of z with characteristic 
length e includes points with the specified conditioning, divided by the scaling factor £2^1/2,0 ^j^^j taken in the limit 
(5 <C e — > 0. These densities are proportional to the correlation function (1), divided by a partition function that 
for a given model is a four-point function of (j)i^2 operators. (The partition function depends on R but not z.) This 
determines the density up to a non-universal factor that depends on the details of how the $1/2,0 operator in (1) 
is regularized. This factor depends on the particular model and is thus not given by CFT. We will ignore it in the 
formulas below. However, in the case of a ratio such as p (see section III) the nonuniversal factors divide out, so the 
conformal result is complete. 



Crossing Probabilities 



We consider the crossing probabilities for the 0(??) model in a rectangle with boundary arcs attached to the corners. 
Inside the rectangle, the boundary arcs join in one of two ways: they may connect the two top and two bottom corners, 
which we call configuration H, or they may connect the two left and two right corners, which we call configuration 
V as shown in figure 1. We use Hh (resp. Hy) to denote the weight of the corresponding configurations for the 0{n) 
model, specifically when the boundary arcs have weight 1. 



H 




V 



FIG. 1: Configurations for the generalization of Cardy's formula. 



We now consider the dense 0(n) model phase. Recall that the Potts model can be exactly mapped onto a represen- 
tation known as the Fortuin-Kastelyn (FK) random cluster model, which replaces the sum over spin configurations 
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with a sum over nearest neighbor bond configurations [16]. This is the basis of the Swendsen-Wang algorithm fre- 
quently used to simulate the Potts model [17]. At the critical point, replacing the sum over bond configurations in 
the FK model partition function with a sum over loops following the cluster boundaries gives a mapping to the dense 
phase of the 0(n) model, where each loop contributes a weight n = ^JQ to the partition function. 

In the FK representation, the <^\^i operators that mark the ends of SLE traces are equivalent to changes between 
free and wired boundary conditions. On a 'free' boundary there is no conditioning and the bond configurations are 
summed over freely as in the bulk, while on a 'wired' boundary all bonds are occupied, so that the boundary sites all 
belong to the same FK cluster. Thus there are two possible interpretations for the correlation function (51): either 
the pairs of operators (?!)i,2(0)(/ii,2(i) and 4'i,2{R)4>i,2{R + i) represent the left and right wired sides of the rectangle 
respectively, or ^i,2(0)(/>i,2(^) and 0i,2(i)<^i,2(-R + i) represent wired bottom and top sides. We will adopt the former 
interpretation (unless explicitly stated otherwise), so the SLE curves follow the boundaries of the FK clusters anchored 
to the left and/or right sides. 

For 0(n) model configurations of type V, the boundary arcs close into two loops that surround the clusters connected 
the left and right sides. Thus the contribution of each such configuration to the FK partition function gains a factor 
and one has 

Zv ^n^Uv . (23) 

For configurations of type H, the boundary arcs close into a single loop around the crossing cluster. Thus the 
contribution to the FK partition function gains a single factor of n and 

ZH = nUH. (24) 

The total FK partition function is therefore 

Z ^ uUh + n^Ilv ■ (25) 

The corresponding result when the top and bottom sides are wired follows on interchanging H and V. 

At this point wc notice a subtlety that will become important when we consider spin cluster densities, namely 
whether the fixed spins on the right side match the fixed spins on the left side. We next construct an FK model 
partition function Z that only samples configurations where the left and right sides are mutually fixed to the same 
spin. This is achieved by adding, to all configurations, a single occupied bond between the left and right hand wired 
sides. Thus, the left and right sides belong to the same cluster, so they are always wired to the same spin. Now for 
the 0(n) model, the boundary arcs must follow the new bond between the left and right edges, so the boundary arcs 
close along the top and bottom sides. The consequences for the FK partition function are as follows. For type V 
configurations, the boundary arcs close into a single loop. Thus 

Zv = nny, (26) 

where Z indicates an FK partition function with the extra occupied bond included. For type H configurations, the 
left and right side are already connected, so the boundary arcs close into two loops. They can be envisaged by placing 
the new bond above the top edge. Then one extra loop follows the top edge of the horizontal crossing cluster but is 
below the new bond, and the other follows the bottom edge of the horizontal crossing cluster, and goes around both 
left and right sides and above the extra bond. Thus 

Zh = Hh . (27) 

The net result for the extended FK partition function is therefore 

Z = n^nH + nnv ■ (28) 

Note that the forms Z and Z are symmetric under exchanges of the labels H and V. This is a consequence of the 
duality transformation that maps the rectangle with left and right sides wired and an extra bond onto a rectangle 
with top and bottom wired and a similar extra bond that is dual to the original bond. 

Now the rectangle with an extra bond is not physical. Thus we need to adjust Z. The necessary factor is determined 
from the explicit form of the critical FK partition function, where each occupied bond contributes a factor of rt. Adding 
an extra occupied bond multiplies each term in the partition function by an extra factor of n. Dividing Z by n, we 
find the FK partition function for mutually fixed left and right sides: 



Z = uHh + Uv 



(29) 



7 



As above, the corresponding result when the top and bottom sides arc mutually wired follows on interchanging H 
and V. 

The result (29) is what we expect. Comparing (25) and (29) we see that the weight of an H configuration is 
insensitive to mutual wiring, while the V weights loose a factor oi Q = n^. 

We now consider the dilute phase of the 0(n) model. This can be related to spin clusters in the Q-state Potts 
model. In contrast to the FK representation, this association docs not follow from direct manipulation of the Potts 
and 0{n) models. Instead one invokes a result known as SLE duality [18], which states that for k > 4 the outer hull 
of an SLE with k has the same fractal dimension as an SLE with k' = 16/ k = (4/7r) arccos(— -y/Q/^), and furthermore 
that these related values of k correspond to CFTs with the same central charge. If the spin clusters can be described 
by an SLE process then it must have parameter k' . 

The boundaries of spin clusters separate regions with a single spin value, say s = 1, from neighboring regions with 
s 1. An SLE that describes these domain walls should separate boundary regions with s = 1 and s 7^ 1 too. For 
the Ising model, the k = 3 SLE trace represents a change between + and — boundary conditions [19]. For the 3-state 
Potts model it has been convincingly argued that the change between s = 1 and s = 2 or 3 (boundary spins freely 
summed over 2 and 3) corresponds to an SLE with k = 10/3 [20]. 

In section V we simulate these dilute phases in a rectangle with s = 1 boundary conditions on the left and right 
sides, freely summing over s 7^ 1 spins on the top and bottom edges. Unlike the FK cluster model, there is only 
one way to fix the boundaries, by fixing both edges to s = 1. This can be implemented by adding an extra nearest 
neighbor bond between the left and right sides exactly as with the FK clusters; the factors found for the mutually 
fixed FK partition function carry through in exactly the same form. Because there is no exact mapping from the spin 
cluster representation to the dilute 0{n) model we cannot fix their relative normalization. Instead, we adopt the Z 
normalization for the spin clusters. 

From the partition functions (25) and (29), we can calculate the probability Ph (resp. Ph) of a horizontal crossing 
FK (resp. either FK or spin) cluster by considering a rectangle with independently (resp. mutually) wired left and 
right sides. This is simply the ratio of the weight of H configurations with the chosen boundary conditions to the 
appropriate partition function. From the above, the weight of H configurations is uIIh for independently or mutually 
wired sides, and the weight of V configurations is n^IIy (resp. Hv) for independently (mutually) wired sides. Thus 

Iljj - nTl^ 

^ UH + nlly' ^" = nun + Uy ' ^^^^ 

For explicit formulas for Il^f and Ily (which make use of results from [21]) see (54), (55), (60) and (61). A result for 
crossing on spin clusters for the Ising model on a circle, which is conformally equivalent to Ph, was obtained in [22]. 
Note that there are other generalizations of Cardy's formula using different boundary conditions (see [8, 18]). 

For n = 1, which is realized by percolation (dense phase) or Ising spins (dilute phase), the distinction between 
mutually wired and independently wired boundary conditions vanishes. For percolation, the partition function for 
either type of boundary condition is constant (i.e. independent of the rectangle aspect ratio) and is renormalized 
to 1 in the continuum limit. Thus, Hh and Uy become Cardy's formula [3] for horizontal and vertical crossings 
respectively, with the property Hh + Hy = 1. 



2. Cluster densities 



The correlation function (1) involves a density operator at the point z, and has contributions from six types 
of configurations {A, B(, Br} and {A, Bf,, Bt} that depend on the location of z relative the the boundary arcs, as 
illustrated in figure 2. We now apply our results to find formulas for the cluster densities in the FK and spin cluster 
Potts models in terms of the weights Hh, Lly, 11^, 11^^ and 11^^. These weights are defined for the 0(n) model, 
specifically when the boundary arcs have weight 1, as for Hh and Ily. In the FK cluster case we find the density of 
sites that belong to the same FK bond cluster as the specified wired edge. In the spin cluster case we find the density 
of s = 1 sites that belong to clusters attached to the specified s = 1 spin edges. 

We now calculate the density Pe{z) (resp. Pr{z), Per{z)) of clusters anchored to the left (resp. right, left and right) 
side(s) of a rectangle with either independently or mutually wired left and right sides. This is equal to some linear 
combination of the weights {n^,!!^^} (resp. {II^jIIb^}, {11^}) divided by the partition function (25) or (29). The 
coefficients of this linear combination are determined by the boundary conditions in the same way as above. The 
boundary arcs of A (resp. Bg, Br) connect as in H (resp. V), so the boundary conditions introduce the same factors 
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PmI th^A 



A 



Bb 




FIG. 2: Schematic drawings of the six configurations associated with (1). We show the boundary arcs and the point z, 
suppress bulk loops for clarity. Note that z must be adjacent to the appropriate side(s) of the rectangle. 



but 



for as they did for Hh (rcsp. the same for Hb,, Hb^, as for Hy). Therefore 



Pr{z) = 
Plr{z) = 



Ha 


+ nlfsj 




+ nWv 








+ nWv 




Ha 







Pdz) = 

Pr{z) = 
Pir{z) = 



(31) 
(32) 
(33) 



As before, the P expressions apply to the independently wired FK configurations, and the P expressions apply to 
both the mutually wired FK configurations and the spin cluster configurations. 

We could also assume that the top and bottom sides are wired, which would result in the same expressions with 
{A, £, r, V, H} replaced with {A, b, t, H, V}, reflecting the (^, ipm) ^ {i/j, ^, 1 — m) symmetry noted in section II A. 

In section V we compare the densities calculated in this section with simulations that measure Pe{z), Pr{z), Pir{z) 
for Q = 1, 2 and 3-state FK clusters and Q = 2, 3-state spin clusters. 

In the next section, we identify the conformal blocks of the six-point function (1) that contribute to the weights 
n^, IIb^ and 11^^ (see (43)-(45)). For exphcit formulas for Uh and Ily see (54), (55), (60) and (61). 



C. Identifying cluster configurations 

From the solutions to the PDEs that govern the correlation function (1) and the prefactor / (10), to go further we 
must identify how the weights 11^, 11^^ and 11^^ contribute to the various conformal blocks G. This can be done in 
an elegant way using Coulomb gas vertex operators, which we will now briefly review; for a more complete treatment 
see, for example, [23, 24]. 

The Coulomb gas representation makes use of a chiral bosonic variable ip{z, z) with action S ~ So + Sc + Sd'- 

So = £^ J {Vffd^x, Sc =^-^ J Uipd^x, Sd ^ a J cos2ipd^x . 

Starting from the free boson action So , adding the complex term Sc that couples the field to the scalar curvature TZ 
modifies the stress tensor. This reduces the central charge of the theory from c=ltoc=l — 24ao^. The field takes 
discrete values, ip G ttZ almost everywhere, due to the term So- This suppresses any rescaling of the parameter g. At 
short distances, the boundaries between these discrete values naturally form a loop ensemble, but at long distances, 
the behavior of the system should still be dominated by So- Reconciling these two scales is only possible if So is a 
marginal perturbation, and this condition lets us relate g and ao, leaving a single free parameter. Thus 

„l-W=<i^^!M, ,34, 

2k 

gives g = 4:/k and 2ao = (1 — 9)9~^^'^- 

Ignoring a subtlety of the zero mode, which is unimportant here, we decompose the boson into holomorphic and 
antiholomorphic components ip[z,z) = '^{z) -\-Tp{z) in order to write chiral vertex operators, 

Vc,(z) = e'^"'^(^) . 
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Chiral operators are sufficient because our problem includes a boundary. The energy flux leaving the system is zero, 
which means (p{z) and 'ip{z) are not independent on the boundary and, by analytic continuation, everywhere else. 
Specifically, in the upper half plane Tp is the analytic continuation of (p from the lower half plane [25] . 

The parameter a is called the charge of the vertex operator, by analogy with two dimensional electrostatics. Non- 
trivial correlation functions must obey a charge neutrality condition. Since Sc effectively adds a non-local background 
charge —2ao to the system the total vertex charges must satisfy J^i — 2ao- The form of charge neutral correlation 
functions is still determined by the original gaussian action So^ 



{JlV^A-^r)) - Z,) 



3>i 



and the charge neutrality condition is entirely responsible for the differences between the Coulomb gas and free boson 
correlation functions. 

The form of the two point function then implies that a given scaling operator has two possible vertex representations 
with charges, a or 2ao — a, related to the conformal weight by 

h = a{a — 2ao) . 

The condition of charge neutrality severely limits the set of computable correlation functions. However, screening 
operators allow one to extend this set of correlations by adding charge without changing the conformal properties 
of the system. Screening operators are non-local with zero weight, and arc formed by integrating weight one vertex 
operators around the singular points of the correlator in a non-trivial way. There are two types of screening operators: 

Q±^ <j> dzVa±(z)- 

The charges and a_ are the positive and negative solutions of 1 = a{a — 2q;o) respectively. Using (34) these are 

a+ = and «_ = --^ . (35) 



We parameterize all other charges as 



± l±r lis 



When r, s S Z+ the charge corresponds to a Kac operator (p^^s- We can generalize to r, s e Z/2, so that the density 
operator has fixed r and s independent of the model parameter k. 

Now we describe how the the vertex operator method helps us to identify the configurations contributing to various 
conformal blocks. The charge neutrality condition for a correlation function with two chiral 01/2,0 ^'Hd four 2 
operators allows a unique set of charges and screening operators: one a^/20' '-"^^ '^1/20' ^'^^''^ ^12^ ^'^'^ Q-S. 
Given (21) we can calculate functions such as (12-16) with the method used in [14], which exploits the common y 
dependence (assumed in [14]) by moving the density operator to the lower boundary, and replacing it with its limiting 
boundary operator. Thus the bulk $1/2,0 operator becomes a boundary magnetization operator (see figure 3) and we 
are left with a five point boundary correlation function made up of one a^3, four a^^^ ^-^^id one 




FIG. 3: The fusion of a bulk magnetization operator to a free boundary has a unique vertex/screening operator representation. 



We choose integration paths for the Q- that entwine neighboring pairs of boundary operators because, with an 
appropriately chosen phase, they yield real solutions for arbitrary values of aspect ratio and model parameter. 

In general we integrate along non-contractible closed paths such as the one on the left hand side of figure 4C. These 
paths can be deformed to run along the real axis except for small circles surrounding the operators. If 4 < k the 
contribution to the integral from the small circles vanishes with the circle's radius, so we can replace the path with 
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^1,2 Cil,2 



B) 



a/,2 



"1.2 W 














a- 



FIG. 4: The boundary five point function requires one a_ screening charge. Integrating the screening operator along the 
boundary as shown gives the five solutions (12-16). 



its component along the real axis. Choosing the paths as in figure 4A reproduces the indicated conformal blocks from 
(12-16). The integration is defined in the upper half-plane, and must be transformed into the rectangle as shown in 
figure 4B. 

For a given correlation function, a conformal block is determined in part by specifying an order in which to fuse the 
operators, retaining a single conformal family from the OPE under each fusion. Mathematically, conformal blocks are 
multiple Frobenius series in small parameters implied by the fusion order with the leading exponents determined by 
the specified conformal families. In the vertex operator formalism pairs of vertices that are not entwined with distant 
vertices have a definite charge when fused together. Such a pair of operators may or may not be entwined by its own 
screening operator circulating around just the pair. In either case, such pairs fuse to a single conformal family, which 
is easily identified by the sum of the fused charges and any entwining screening charges. Thus a given set of screening 
charge integration paths in an arbitrary correlation function is simply identified with a conformal block. 

Next, as was done in [14] for percolation, we use the vertex operator formulation to determine the contribution 
of the conformal blocks to the configurations of interest. First, for the integration paths in figure 4 we deduce the 
fusion structure of each block to determine the contributions from the three configurations illustrated in the first row 
of figure 2. Recall that these weights are defined for the 0(n) models, with the weight of the boundary arcs being 1. 
These conformal blocks can be identified with sets of boundary conditions similar to those in section II B, except that 
we may leave unattached arcs instead of boundary loops. 

First consider the block Gi in the limit where the magnetization operator goes to the lower left hand corner. This 
traps the screening operator as well, for a total charge of Q;j~2 + ctf^ + a_ = af^, which is the charge for a three leg 
operator. An n-lcg operator, associated with c^i^i+n, represents a point where at least n non-contractible boundary 
arcs (arcs that cannot be smoothly deformed to a point because of some constraint on their behavior) each with weight 
1 are attached to the boundary. The only configuration with three or more non-contractable boundary arcs emanating 
from this corner in that limit is Br', therefore this configuration must be the sole contribution and Gi ~ 11^^. (Note 
that this proportionality disregards all y dependence, which is contained entirely in the prefactor /(^, ■)/), m).) 

The block Gu is more subtle. Fusing the right hand corners together yields a charge of 2aj~2 = aj^g, implying a 
two-leg operator; this excludes configurations which have contractable boundary arcs attached to the right side. 

The remaining contributions to Gu are fixed by considering the left hand edge, which has a net charge of 2Q;j~2+Q;_ = 
indicating that the one leg operators at the corners fuse to the identity family. If the arcs attached to this edge are 
non-contractable (as with A configurations) then they must be attached to other 1-leg operators. These operators 
will fix the arc weight at 1, so the total weight does not change. If the attached arc is contractable, as with Br 
configurations, then bringing the operators together would form a loop with weight 1, and a compensating factor of 
n is needed. The net result is that Gu ^ ILa + nTls,.- The same result follows by adding a boundary arc between 
each pair of one-leg operators that fuse to the identity family. 

When we consider the block Gui the top edge has net charge 2Q;j~2 + a- = and fusion is through the identity. 
Since there are no other restrictive fusions this block contains contributions from all three configuration types. The 
top edge of type A configurations has a contractable loop, while Br and Bg have non-contractible loops attached 
to the top edge, thus Gui ^ nllA + 115^ + Ubc- Again this result can be obtained via a boundary arc, this time 
connecting the top corners. 
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The associations for Giv and Gy foUow from mirror symmetry and for Gvi by definition. The complete list is 



Gi ~ Hb,. , (37) 

Gii ~ + uUb^ , (38) 

Giii ~ hRa + Hb,. + Hb, , (39) 

Giv - TIa + tiUb, , (40) 

Gv - Hs, , (41) 

Gvi - . (42) 



The relative normalization in (12-16) emerges naturally in the vertex operator formalism, since the integral expres- 
sions for the different conformal blocks have different integration paths but the same integrand. Hence the proportion- 
alities in equations (37-42) all omit the same factor. This is an advantage of using vertex operator over the differential 
equation analysis of [1], which would require comparing several limiting cases to fix the relative normalizations. 

Using the expression (10) for the common pre-factor /(^, ip, m) we thus find explicit expressions for the weights of 
the three configuration types A, Br and Bg, as 

Ua = /(e,V',™)Gvi(^,m) (43) 
Ub,. = f{C,^,m)Gi{^,m) and (44) 

ns, = /(e,V,™)Gv(^,m) . (45) 

Explicit formulas for / are given in (lO)-(ll); and for Gi, Gy and Gvi in (12), (16) and (18), respectively. Recall 
that 5 is independent of y and the common factor /(^, ip, m) contains (via V') all of the y-dependence in these weights. 
This plays an important role in section III. 

Returning to the functions Gj-Gy, we can also derive relations analogous to the four-point crossing matrix directly 
from the vertex operator formalism. The integral of the screening charge around a small loop in the bulk containing 
no operators is zero. Deforming the contour to run along the boundary gives the blocks as described above, while 
picking up a complex argument of each time the contour passes a boundary charge a. Thus 



= e«'^*/"Gi + e'^'^^/'^Gii + Gm + e.-'^'^^/^Giy -I- e-^'^'/'^Gy . (46) 

Making use of (19) then allows us to write (46) as two real equations 

2Giii (2-n2)Gi+nGii + nGiy + (2-ri2)Gy , and (47) 

nGi + Giy = Gii + nGy . (48) 



which also follow from (37-41). These two equations reflect the fact that the five solutions live in a three dimensional 
solution space, which is consistent with the number of configurations in the top row of figure 2. 

In the foregoing, we have simplified our figures and discussion by considering linear integration paths along the real 
line only. However, our results also hold for the more general paths entwining the operators required for k < 4. Thus 
our solutions and the identifications of their fusion channels extend automatically to the dilute phase of the 0(n) 
model, and equivalcntly the spin cluster behavior of the critical Q-statc Potts models. 

While the Coulomb gas formalism has considerable predictive power, our differential equation analysis has the 
advantage of specifying the form of the prcfactor (10), and determining the dimensionality of the entire solution space 
without making a priori assumptions about the properties of the density operator. In addition, (21) does not follow 
simply in vertex operator formalism. 

The three cluster configurations illustrated in figure 2 account for the three dimensional solution space for 
F{£_,il},m) = G{^,m) discussed above. In addition there are three solutions that occur when the top and bot- 
tom rather than left and right sides are wired. These additional solutions are implied by the {^, ip, m} O — w} 
symmetry (22). This exhausts the set of solutions that are consistent with the conformal weights and null state 
conditions of the correlation function (1). 

III. FACTORIZATION BEHAVIOR 

In this section we consider a universal ratio that demonstrates a factorization of certain higher-order correlation 
functions (or the equivalent densities) in terms of less complicated correlation functions (or densities) . The ratio was 
originally defined at the percolation point (k = 6), and is discussed in [11-14]. Here we generalize it to a set of 
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two-dimensional critical points. The generalized ratio is universal, constant (except near the sides of the rectangle), 
so that factorization occurs, and independent of the vertical coordinate y. 
In [14] the universal quantity 

( m P£r{z,R) 

Pperc ix,R)= , (49) 

^/Pr{z, R)Pe(z, R)Ph[R) 

was considered for critical percolation clusters. Here Ph{R) is the horizontal crossing probability given by Cardy, and 
Pir{z,R) (resp. Pr{z,R), Pe{z,R)), already introduced in subsection IIB2, are the densities of clusters connecting 
to both (resp. right, left) vertical sides of the rectangle. (These densities are given in terms of Ha, and by 
(31)-(33) with n = 1. Note that in this case Hh + Hy ~ 1. Note also that all nonuniversal factors divide out of 
the ratio, so it has the same value if cluster densities in a specific model of percolation are used in place of the CFT 
results.) 

The y-independence of Pporc is discussed in [14] and herein. 

In the limit i? — >■ cx) Pperc is constant and Pir exactly factorizes in terms of two- and three-point functions. We 
showed that in fact /5perc(oo) = 1.0299. . . is the fusion coefficient for three two-leg operators. See [11-13] for more 
details. 

In a finite rectangle PpGrc(2;, R) is not quite constant, though it is independent of the vertical coordinate, as implied 
by (21). Instead ppcrc{x,R) equals 1 at a; = and i?, and decays exponentially, on a distance scale of the height of 
the rectangle, towards Pperc(oo) within the bulk. Thus Pperc is constant to within 3%, and the factorization holds to 
good approximation for arbitrary x and R at the percolation point. 

The results of the previous sections allow us to generalize Pporc to arbitrary k. To begin, we introduce the notation 
le := f/)! 2(0)^1 2(1)' '■= 01 2(^)'/'i + i)j ^-i^ti cr(-Zj^) '■= ^i/2,o(-2^j ^) for the wired left and right sides of the 
rectangle and the magnetization operator at interior point z respectively. We then generalize the definition from [14] 
as 



/ - m / (2^f[i,3]0-(2,2)[i,3]2:r)2(2:f[i4]2:r) 
p{z, z, R) = J — — — — -— — — — — — ; (50) 

the motivation for this particular choice is explained below. The subscripts index the propagating channel between 
sets of operators, thereby identifying the conformal block. The new correlations introduced here specify weights, when 
properly normalized these weights become densities (if z dependent) or probabilities (if not). 

The correlation functions (Iilr) with four (j^l 2 operators, in the cases where they can be understood in terms of 
critical clusters, generalize Cardy's horizontal percolation crossing probability, as discussed in subsection II B. The 
configurations in figure 1 determine the blocks (see (59-60) below), which are given by 

(Z,[^.]Z,) = {c^l2i^)<l>UiO)[x]'l^l2iRWi,2i^ + R))n (51) 

n'^ \w'(z(v ■]m'^^'^ 

= H n ^ IIL {Mn,)Mni)[x]Mn2)M^,))w (52) 

£j^0 (16ei£2£3£4)^'^ 

= {K'f'^^^G^{l~m) . (53) 

The covariance prefactor in (52) is computed in [1]. The conformal blocks for this four point function were given in 
[21] for a multiple SLE process. Due to the factors arising from the corner operators we use a slightly modified form 



Gi,i(m) = 2F1 (2-— ,l--;2- 



K K 



r(12/K-l)r(2-8/^) ("i 4 4 

^'''^'^> - r(8/^)r(i-4/^) ™ 



(54) 
(55) 



We next deduce the form of the conformal blocks for the six-point function using the methods from section II C and 
(43)-(45). Two (f>i^2 operators that fuse in the 0i,3 channel cannot represent the two ends of a single boundary arc. 
Also, while two 0i.2 operators that fuse in the identity channel do not place any conditions on which configurations 
contribute to the block, they do imply an extra factor of n for configurations in which they are connected by a single 
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arc. It follows that 



It ^{ ~ 




[l,3]Ir; - 


lU 




(56) 




.^") 


[l,3]Ir) = 


Ha 


f tiIIb^ = /(^, V', m)Gii(^, m) 


(57) 




.^-) 


[1,1]^'-) = 


Ha 


f nllB^ = /(^, ip, m)Gw{^, m) 


(58) 




(If 


[l,l]Ir) = 




+ nUv =(i^')"''''Gi,i(l-m) 


(59) 




(If 


[l,3]Ir) = 






(60) 


One may also write 
















Hy = 




)«''^'^Gi,3(m) . 


(61) 



The method of determining these contributions via n^s uses 0(7i) loop model concepts, and thus applies when 
8/3 < K < 8. But conformal blocks arc uniquely identified by the exponents in their various functional expansions, 
so once identified the results extend to all k > 0. 
The expression for p{z, z, R) thus becomes 



Gvi(e,m)2Gi,i(l 



Gii(e,m) Giv(e,m) Gi,3(l-m) 



(63) 



Note that the independence of p{x^ R) from the vertical coordinate y is manifest. This feature may be traced back to 
(9). 

The particular choice (50) for p{x, R) ensures that each conformal operator in the numerator has a counterpart 
in the denominator. This guarantees that we may replace each correlation function by its corresponding density, 
since all the nonuniversal scaling factors cancel regardless of the regularization scheme, so that the resulting ratio of 
densities is a universal conformally invariant quantity. We emphasize that this generalization is not unique. It does, 
however, minimize the number of conformal blocks that appear in the definition, and in that sense may be considered 
the natural generalization of Ppcrc- 

The left hand sides of equations (31)-(33) and (30) then allow us to write the universal ratio p in terms of the 
densities Pi and the crossing probability Ph for FK clusters in rectangles with independently wired vertical sides as 

R) = , '^"^''''^ , (64) 

^Pl{z,R)Pr{z,R)PH{R) 

which naturally generalizes ppcrc in (49) to 8/3 < k < 8. For spin or FK clusters in rectangles with mutually wired 
vertical sides one may write p in terms in terms of the densities Pf (z, i?), Priz, i?), Pi,r{z, R) and crossing probability 
Ph{R) by using the right hand sides of (31)-(33) and (30), respectively. By inverting these equations, we find 



, „^ P,MRY[n^-{n^-l)PH{R)] , . 

'^^ ' ^ ^ [{n^^l)Per{z,R)-n^Pi{z,Rmn^^l)hr{z,R)-n^Pr{z,R)]P„{R)- ^ ' 

When Q — 2, n — 1 (dilute) and the expression (65) for p reduces to (64). This occurs because the left and right 
sides are wired + (resp. — ) while the top and bottom sides are wired — (resp. +) when Q = 2. The consequent +/— 
duality eliminates any distinction between closing loops along the left and right sides of TZ versus closing loops across 
the top and bottom. One could define p differently so that (65) was simpler, e.g. by (64), with each P replaced with 
a P. However doing so results in an expression that does not reduce to a single OPE coefficient in the limit 7?. — )■ oo 
(see (67) below). 

The ratio p(x, R) is normalized so that p = 1 when we take x to or R. In terms of critical clusters, this means that 
as z goes to the left hand side it becomes vanishingly likely that z is not connected to the adjacent wired boundary 
and — > 0. In addition, 11^ — > a-nd n^;, — > By, so that indeed p — > 1 as claimed. An analogous argument 
holds as z goes to the right hand side. 

We can also take the long rectangle limit R oo with <C .t ^ i?. Then the existence of a horizontal spanning 
cluster is exponentially unlikely so that B^ ^ B^^, B^^, and Hh ^ By. Furthermore, we can replace the distant left 
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and right sides with point operators, and due to the invariance with respect to the vertical direction we can assume 
the magnetization operator sits on the boundary. Thus 

p(a;,i?->oo) = / TT^n ^ ^^^^ 



^n(0i,3(a;)0i,3(i?))(0i,3(O)0i,3(a;))(</'i,3(O)0i,3(i?)) 



\ -Wr(^)r(i^)r(f)= 



(67) 
(68) 



where the final line uses the known expression for the OPE coefficient as given in [26]. (Recall that n^n) is given in 
(19).) Note that the value of C^'3.1 3 increases monotonically from to 5-\/2 as k decreases from 8 to 2. Some specific 
values are given in the captions to figures 8-11. For percolation (k = 6) Cl'^.^ 3 is close to 1, and the factorization is 

similar anywhere in the rectangle, since p = 1 at the edges. When C^'3.]^ 3 is far from 1, however, the nature of the 
factorization is different near the edges and far from them. 

Note that the long rectangle limit, as described above, fails when k = 8/3. This value of k corresponds to the dilute 
phase with n = and from (62) it follows that p{x, i?) = 1, in apparent contradiction to (68). This is because the two 
limits K — > 8/3 and i? — > 00 do not commute; vertical crossings are suppressed as n — >■ 0, while horizontal crossings 
are suppressed as i? — > 00. This problem does not occur in the dense phase with n = (k = 8). Here the limits do 
commute and p = 1, because the curves are space filling, so there is no suppression of horizontal crossings for large R. 

The form of this ratio for the (5=1,2 and 3 state Potts models (for either FK or spin clusters) is given in figures 6, 
8, 9, 10 and 11 (solid curves) for various values of the aspect ratio R. Note that for i? = 2 and 3 there is virtually no 
difference between the p(a;, R) curves for a given Q value when < x < 1, although the spin clusters deviate slightly 
more. This is consistent with the observation made in [14] for percolation {Q = 1) that the semi-infinite strip result 
p{x, 00) is a good approximation to p{x^ R) for Q < x < R/2 and i? 3> 1. 

Note that taking the limit i? — > 00 with ^ a; ^ i?, in which case p(a;, R) becomes a constant, is equivalent to 
bringing the pair of 4>i^2 operators in (1) on the left hand side of the rectangle together, and likewise the pair on the 
right. Applying this to (X^ [1 3](t(z, z)[i^^Xr) one obtains a four-point correlation function with two 0i^3 operators; the 
other correlation functions in (50) can be treated similarly. The result generalizes our exact percolation factorization 
formula for the densities of clusters anchored to two points [12] to all the critical models mentioned here. By keeping 
X near 0, one could similarly obtain results for the densities of clusters anchored to one interval and one point. We 
will examine the consequences of these observations elsewhere. 



IV. DENSITY OF PERCOLATION CROSSING CLUSTERS 



We now restrict our attention to k = 6, or critical percolation. We return to percolation, not just for its intrinsic 
interest, but because we are now in a position to exploit the property known as locality in the SLE literature [8]. This 
property allows us to use the results from the subsection II C, with wired boundary conditions on the vertical sides, 
to derive an explicit expression for the density of clusters that cross a rectangle with free boundary conditions. 

With respect to a growing SLE hull, 'locality' means that as the curve grows in time, its statistics are completely 
independent of boundaries up until the hull actually touches one. Percolation exhibits this property because the 
decision to include or exclude a link is determined by a purely local random variable and not by the state of neighboring 
sites. This implies that the bulk sites all contribute a fixed amount to the overall weight of a configuration regardless 
of the state of the boundary sites. 

It is this insensitivity to the boundary conditions that allows us to generalize the idea of wiring. Suppose we 
condition boundary sites in some way, and find the density of clusters attached to certain segments, as above. If we 
change the conditioning and measure the density again, then the two densities are directly comparable due to the 
fixed contribution of the bulk sites to the overall weight. 

This is intrinsically different from a system without locality. In general conditioning a set of boundary sites changes 
the degrees of freedom of the attached clusters, so boundary conditions are a defining characteristic of the system. For 
a system with locality we needn't think of boundary conditions as fixed parameters; instead they provide a variable 
means of probing the connectivities of clusters within the system. 

We begin with the three configurations illustrated in figure 2, and proceed to decompose the contributions to Ilyi, 
as in figure 5. We distinguish configurations according to whether the point z actually belongs to a crossing cluster, 
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FIG. 5: Percolation configurations discussed in tiiis section. We use fixed boundary conditions on botfi (resp. riglit, left) vertical 
sides to calculate the densities corresponding to Br, Be, A (resp. Cr, Ce)- Due to locality, however, this decomposition of the 
densities holds regardless of boundary conditions. The two top rows represent the disjoint subsets that span the configuration 
space. The bottom row represents quantities that we calculate here which are composites, as indicated by the paths from their 
constituents in the upper rows. 



as in configuration type , or whether it connects to the crossing cluster only through the wired boundaries on the 
right or left hand side, as in type Ar and Ai respectively. Our goal is to determine the density for configurations of 
type Ax, since these actually connect both vertical sides of the rectangle, even when they are not wired. They also 
constitute the horizontal crossing clusters for a rectangle with free boundary conditions on all sides. Thus 



Pa ^Pa 



Pa, + Pa 



(69) 



For percolation, the expressions for the conformal blocks G in (12)-(18) are equivalent to those calculated in [14], 
however (10) gives us an expression for the prefactor /(^,'0,to), which contains the explicit vertical dependence not 
previously available, thus allowing a complete calculation of the densities. 

Now we consider the correlation functions given by wiring either the left or right hand sides of the rectangle. In 
[11] we calculated the density of percolation clusters at a point w = u + iv vn the upper half-plane that are attached 
to an interval on the real line / = (ui, U2), 



Piiw) 



r5/4M2Im 



1/4 



1/3 



1 



(w - Ui){w - U2) 



(70) 



(w - Ui){w - U2) 

1, so the weights found in [11] and also those from 



Note that for percolation the normalization factor Hh + nliv 
(1) are in fact densities. 

We transform Pi{w) into the rectangle TZ using the mapping w{z) = to sn {zK'\m) from [1]. The interval {—00, 0} 
will map to the left side and {m, 1} maps to the right side, thus the corresponding cross-ratios are 



Vr 



— = exp (4i arg(sn (z^^T'Ito))) 



{w — m){w — 1) 



= exp (4i arg(de {zK'\m))) 



(w — m)(w — 1) 

This gives, if the right or left side of the rectangle is wired, respectively 



(71) 
(72) 



PcAz,z) = 2^/%K')^/^^sm^/^[siTg{sn{zK'\m))] 



Im 



sn {zK'\m) 



|sn {zK'\m) cn {zK'\m) dn (z/^T'Ito)! 



Im 



sn {zK'\mY 



-5/48 



-5/48 



|sn {zK'\m) cn {zK'\m) dn {zK'\m)\ 



(73) 



(74) 



The coefficient 2^/^ = 2'*^ indicates the presence of a corner operator, as in (11). 
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Using identities for elliptic functions and (10) we can rewrite these quantities in terms of the real coordinates (8): 



PcAz,z) = f{^,^p,m) [dc^{xK'\m) + mcd^{xK'\m) - dn^{yK'\l - m) - mnd^{yK'\l - m)] 



1/6 



2l/3(-^'-)5/48 



(1-TO^2)2 



(1 - (1 - m)tp'^)' 



e(l- 0(1-0 i>il - - {I - m)^P) 
(1 - (1 -?n)V')2 +m 



-4 



-11/96 



1 - 



1 — (1 — m)?/' 



1/6 



Pcf {z, z) = f{£_,ijj,m) \jis^{xK'\m) + msT?{xK'\m) — dn^(?;_ftr'|l — to) — mnd?' {yK'\l — to)] 



1/6 



2l/3(-^')5/48 



(1 - TO^) 



2\2 



(1 - (1 - to)V'^) 



2\2 



e(l-e)(l-mC) 7/;(l-V)(l-(l-™)V') 



-11/96 



TO^ 



(1 — (1 — m)?/>)^ + TO 



1 — (1 — m)i/; 



1/6 



Note that on use of (22) and the invariances of f{£^,ip,m) [1], the (final) expressions for Pp,. and transform 
properly under mirror symmetry about x — R/2 or y = 1/2. 

With these expressions for the densities Pc,. and we can isolate the density of crossing clusters Pa^ . Now 



Pc^ = Pb^ + Pa,. + Pa^ 
Pc, = Pb, + Pa, + Pa^ 



and 



(75) 
(76) 



Note that in (75) and (76) the left and and right hand sides involve quantities originally defined with different boundary 
conditions. Only when there is locality does this difference become irrelevant. 

We now check that the normalizations of the various quantities are consistent. Using (58), (76) and figure 5, we 
find /Giv — Pc'e = Pa^- The latter is the density of points z that connect to the right hand side but not to the left, 
in configurations with a crossing cluster. As x — )■ 0, our expression for this vanishes, as it must. Using (57) and (75) 
we come to the same conclusion for / Gu — Pc,. = Pai as a; — )- i?. 

Thus we arrive at an explicit formula for the density of clusters that cross the rectangle horizontally regardless of 
wiring at the sides (equivalently, cross a rectangle horizontally with open boundaries). 



Pa^ - Pc. + Pc, -Pa- Pb^ - Pb, 
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1 — (1 — m)^ 



(77) 



1/6 



(78) 



(1-m) 



m-m-mO]'^' 'V3'3' 3' 3 



m, 



(The invariance of the last term in the bracket in (78) under mirror symmetry about x — R/2 follows from (22) and 
standard results for the Appell functions.) 



V. SIMULATION RESULTS FOR p{x,R) AND Pa^ 



In this section, we present numerical results that verify the predictions (43-45) via calculation of the universal 
ratio p{x, R) (see (62)). Here we simulated the Q — 1, 2, and 3-state Potts models on a square lattice in rectangles TZ 
of aspect ratios R — 1/2, 1, 2, and 3 using the Swendsen Wang (SW) algorithm [17] at the critical bond activation 
probability Pc{Q) = y/Q/{l + \/Q). Both FK and spin clusters were considered, and we found good agreement 
between the predictions of this paper and simulation. We also verify (78), the density Pa^ of percolation clusters 
that cross horizontally with free boundary conditions via simulations. 
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FIG. 6: Upper graph: Numerical results for p{x,R) for percolation {k = 6) against the CFT prediction (solid curve). The 
data points are simulation values of p(x, y, R) averaged over the ^/-coordinate with R fixed. Simulations were performed on a 
rectangle of about 10^ square lattice spacings. Notice that both the theory and data plateau at C^'f-i 3 = 1.02993 (dashed line) 
given in (68). The lower graph shows, at fixed x,R, the standard deviation of measured p{x,y,R) values from their average 
along the height of the rectangle. 
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FIG. 7: Standard deviation of measured p values for figures 6 and 8-12 at fixed x and R from the measured averages, as in 
the lower graph of figure 6 (percolation). The average of these standard deviations over x is shown. In each case, the graph of 
standard deviation vs. x resembles the the lower graph in figure 6, but the magnitude varies. 
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FIG. 8: Same as upper figure 6, except for Q — 2 FK cluster densities = 16/3). Here the plateau is at CYt \ 3/2^''* ~ 1.07871 
(dashed line), from (68). 
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A. Q = 1: Percolation 

1. The universal ratio p(x,R) 

When Q = 1, the rectangle TZ is filled by a single spin cluster in one spin state. The FK clusters within it are bond 
percolation clusters, and because Pc(l) = 1/2, the SW algorithm is equivalent to simulating critical bond percolation. 

The boundary conditions are either "fixed" ("wired") or free ("open"). All bonds on a "fixed" interval are activated; 
on a "free" interval they are activated with probability 1/2 (just as in the interior of 7?,). The insertion of (j)i^2 changes 
the boundary condition from fixed to free, so the boundary conditions on the sides of TZ alternate between fixed and 
free. We choose to wire the left and right sides, so the top and bottom are free. 

In our simulations, we measured the density Pi{z,R) (rcsp. Pr{z, R), Pir{z, R)) of clusters touching z G TZ that 
are anchored to the left (resp. right, left and right) side(s) of TZ. The density was found numerically by dividing the 
number of samples where z is connected to the left (resp. right, left and right) side(s) through activated bonds by the 
number of samples. We also measured the probability of horizontal cluster crossings Ph{R)- The length and width 
of TZ were chosen so that the area of its interior is approximately equal to 10® square lattice spacings. When i? = 1/2 
or 1, 5 X 10® samples were generated, and when i? = 2 or 3, 15 x 10® samples were generated. 

In figure 6, we plot simulation results for ppcrc{x,R) and the CFT prediction in (49). We expect Pperc to be 
independent of its y-coordinate, so in the first plot, we average the measured Pperc values over the y-coordinate for 
fixed X and plot the average as a function of the x-coordinate, as well as its standard deviation from the measured 
values. When x nears oi R, p decays exponentially to 1. 

Some of our percolation results overlap with the simulations in [14]. However, we include them here for ease of 
comparison with our results for Potts models at Q = 2 and 3, which are new. 

B. Q = 2 and 3-state Potts Models 

Here, two different types of clusters may be considered: FK clusters and spin clusters. In the continuum limit, 
the boundaries of FK clusters anchored to the sides of TZ are SLE curves in the dense phase with speed k related to 
Q — r? via (19). When Q = 2, k ^ 16/3, and when Q = 3, k = 24/5. On the other hand, the boundaries of spin 
clusters are SLE curves in the dilute phase with speed k' . As described in subsection II B, k' = 16/k, where k, is 
the speed of the SLE curve for the corresponding FK cluster boundaries. When Q = 2, k' = 3, and when Q = 3, 
k' = 10/3. 

The boundary conditions to be used will depend on which type of cluster, FK or spin, we wish to study. This is 
because the boundary condition change induced by the insertion of c/)i_2 operators at the corners of the rectangles is 
different in the dense phase from that in the dilute phase. Consequently there are two different ways to implement 
the boundary conditions, as discussed below. 

1. FK clusters: the universal ratio p{x,R) 

In simulations of the Q-state Potts model in a rectangle TZ with independently wired left/right sides, we measured 
the density Pe{z, R) (resp. Pr{z, R), Pe,r{z, R)) of FK clusters at lattice site z = x + iy G TZ connected to the left (resp. 
right, left and right) side(s) of TZ. These densities were found by dividing the number of samples where the lattice 
site z touches an FK cluster anchored to the left (resp. right, left and right) sides by the number of samples. In the 
same simulations, we also measured the probability Ph{R) of a horizontal FK cluster crossing (which can occur only 
in samples where the left and right sides have the same spin), given by (30). The densities are given in (31)-(33), and 
the probability Ph{R) in (64). 

As discussed in subsection II B, 0i^2 changes the boundary condition at the corners of TZ from fixed to free in the 
dense phase. In the FK representation, "fixed" implies that all bonds along the left and right sides of TZ are activated 
in each sample, and "free" implies that bonds along the top and bottom sides are activated with probability Pc{Q). 
To ensure that these boundary conditions are respected, we modify the Swendsen-Wang algorithm [17] by laying a 
bond with probability 1 between all adjacent sites of the left and right sides before each update. 

The system size and sample size were chosen differently for Q = 2 and Q = 3. When Q = 2, the length and width 
of TZ were chosen so that the area of its interior is approximately equal to 10® square lattice spacings. 5 x 10® samples 
were generated when R = 1/2 or 1, and 15 x 10® samples were generated when i? = 2 or 3. When Q = 3, we noticed 
that the results for the universal ratio p (64) in a rectangle of 10® square lattice spacings were noticeably lower than 
the theory prediction, so in this case we quadrupled the area of TZ to about 4 x 10® square lattice spacings when 
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i? = 1, 2 or 3, which improved the agreement. For Q = 3, 5 x 10^ samples were generated when R = 1/2, and 15 x 10^ 
samples were generated when i? = 1, 2 or 3. 

In figure 8, we plot values of p measured in our simulations for (5 = 2 (k = 16/3) and observe good agreement. We 
expect p to be independent of its y-coordinate, so we average p over the y-coordinate for fixed x and plot the average 
as a function of the x-coordinate. In figure 8 (resp. 9), we plot values of p measured in our Q = 2 (k = 16/3) (resp. 
Q = 3 (k = 24/5)) simulations. Also for fixed x, we measured the standard deviation of the measured p{x,y) from 
the corresponding measured average over y. The average of these standard deviations over all x values is shown in 
the table in figure 7. For small R, these deviations are less than 0.013%, and for larger R, they arc less than 0.63%. 
Notice that when R is large (w 3), p is roughly constant except near the sides of the rectangle and equal to the fusion 
coefficient divided by ^/n (68). When x nears or R, p decays exponentially to 1. 

2. Spin Clusters: the universal ratio p{x,R) 

In simulations of the Q-state Potts model in TZ with mutually fixed left/right sides, we measured the density Pe{z, R) 
(resp. Pr{z, R), Pe.r{z, R)) of spin clusters at lattice sites z = x + iy G TZ connected to the left (resp. right, left and 
right) side(s) of TZ. These densities were found by dividing the number of samples where the lattice site z is in a 
spin cluster anchored to the left (resp. right, left and right) side(s) by the number of samples. We also measured the 
probability of horizontal spin cluster crossings Ph (R) ■ 

As mentioned in subsection II B, 0i^2 in the dilute phase changes the boundary condition at the corners of TZ from 
spin type A (left and right sides) to an unbiased mix of spins that are not A (top and bottom sides). Notice that 
both the left and right sides are necessarily mutually wired, in contrast with FK clusters. 

Because we cannot allow a spin cluster of type A to touch the top and bottom of the rectangle, how to maintain 
the boundary condition change at the corners with each update is different from that of FK clusters. Suppose that 
we start with some arbitrary spin configuration inside TZ that satisfies the necessary boundary conditions. To update 
the system, wc activate all bonds on the left side and grow the FK cluster anchored to it via the SW algorithm. Note 
that this FK cluster cannot touch the top and bottom since these sides do not contain type A spins. If the FK cluster 
strikes the right side, the whole right side is absorbed into it because the right side is wired. If it doesn't, then we 
repeat this process on the right side, and we update the spins in the FK clusters anchored to either side to type A. 
Next, we grow the FK clusters anchored to the top and bottom, and their spins are updated to a type other than 
A with uniform probability. Finally, we grow an FK cluster from any site in the interior of TZ whose spin is not yet 
updated, with no restriction on the updated spin type. 

For given Q and i?, the rectangle size and number of samples were chosen in the same way as for Q = 2, 3 FK 
clusters. When Q = 3, we noted that the results for the universal ratio p (65) in a rectangle of 10^ lattice spacings 
were noticeably greater than the prediction, in contrast to FK clusters. So to obtain better results for Q = 3, we 
again quadrupled the area of TZ when i? = 1, 2 or 3. 

In figure 10 (resp. 11), we plot values of p measured in our Q = 2 (k = 3) (resp. Q = 3 {k = 10/3)) simulations. 
Notice that the data at i? = 1/2 is too small to see. It is plotted on a more visible scale in (12). As done earlier, we 
plot averages of the measured p over the y-coordinate for fixed x. Standard deviations of the measured p{x, y) from 
the corresponding measured average over y are shown in the table in table 7. For small _R, these deviations are less 
than 0.011%, and for larger R, they are less than 15%. We comment on possible reasons for large deviations in the 
latter case below. Again, when R is large (~ 3), p is roughly constant and equal to the fusion coefficient divided by 
(68). When x nears or R, p decays exponentially to 1. 

The standard deviations for spin clusters at large R shown in table 7 are noticeably greater than their FK coun- 
terparts. Indeed, figures 10 and 11 show deviations from the mean value of p{R = 3) greater than those in the 
corresponding FK cases. This may be attributed to the increase of conformal weight hi^^. In the large R limit (see 
(68)), p is given by two- and three-point correlation functions of 0i_3 operators. For spin (resp. FK) clusters at Q = 3, 
hi^3 = 7/5 (resp. 2/3), while for Q ~ 2, /11.3 = 5/3 (resp. 1/2). Thus the correlation functions decay more rapidly 
with distance in the spin case, making the events Pi, P^, and Pi ^ less likely than the events P^, P^, and P^^j. in the 
FK case, so that simulation errors increase. Assuming that the ratio of the deviation for the spin case to the FK case 
is roughly equal to corresponding ratio of p values, we find (for x near i?/2) a ratio of order 10 in either case which 
is consistent with our results. 

Figures (11) and (12) shows that the Q = 3 case (k = 10/3) also exhibits these features. The most noticeable 
difference between the Q — 2 and Q = 3 cases is the considerably larger error in the latter when i? = 3. Because 
the system size of the latter is quadruple that of the former, the (5 = 3 run time for a single sample is quadruple 
the Q = 2 single sample run time. Also because more samples need to be generated in order to achieve a reasonable 
accuracy, we were unable to achieve the same degree of accuracy as with the Q = 2 case. The (5 = 3 data fluctuates 
about the theory curve, so the large error is most likely attributable to noise that can be reduced with more samples. 
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C. Density of percolation crossing clusters with free boundary conditions 

The density Pa^ of horizontal percolation crossing clusters with free boundary conditions on all sides of TZ is given 
by (78). The data (resp. prediction) was normalized by dividing by the measured (resp. predicted) value in the center 
of the rectangle. In order to simulate Pa^ for percolation, we carried out a simulation similar to what was done in 
[14], in which we "grew" clusters from all the sites at one edge of the rectangle. For those clusters that reached the 
opposite edge, we added all wetted sites to an array that tallied the number of times each site of the lattice was 
visited. Dividing by the total number of runs gave the density. If one cluster crossed, we did not include the sites 
of any other, non-crossing clusters, as would be the case for fixed ("wired") boundary conditions. We considered 
square and rectangular lattices, results for 128 x 128 and 128 x 256 are shown. For larger aspect ratios, the statistical 
significance dropped because of the lower probability of finding a crossing cluster. For each, we carried out about 10^ 
samples. 

In figure 13 we show a comparison between the theory and the measured values. The contour lines arc at heights 
0.5, 0.6, 0.7, 0.8, 0.87, 0.93, 0.96, 0.98, 0.994, 0.999, from the outside inwards. The deviations are quite smaU. For 
instance, for i? = 2 and various fixed x values, we found the following relative errors (absolute value of the difference 
between predicted and measured value divided by the predicted value) averaged along the y coordinate: (x, error) = 
(1, 0.0094), (0.5, 0.0057), (0.25, 0.0038), (0.125, 0.0036). And for various fixed y values: {y, error) = (0.5, 0.011), 
(0.25, 0.0055), (0.125, 0.0030), (0.0625, 0.0021). 

VI. SUMMARY 

In this section we summarize the main results in some detail, and reproduce the main formulas. The research done 
here exploits, for a variety of two-dimensional critical models, the solutions for a chiral six-point correlation function 
in order to specify the density of critical clusters anchored to one or both vertical sides of a rectangle with wired 
boundary conditions on those sides and free boundary conditions on the horizontal sides. These models include the 
critical 0(n) loop models in both high and low density phases and both Fortuin-Kasteleyn (FK) and spin clusters in 
the critical Q-state Potts models. We also determine the density of percolation crossing clusters (from one side of a 
rectangle to the opposite) with free boundary conditions. 

Subsection II A reviews the solutions found in [1] for the correlation function (1) 

C{Z) = ((/)J,2(O)0^,2(i)*l/2,o(^,^)'/'t2(^)'/'t2(^; + i)>7^ , 

in the rectangular geometry TZ -.^ {z = x + iy £ C\0 < x < R,0 < y < 1}. Here z ^ x + iy, C also depends on 
the parameter m that specifies the aspect ratio of the rectangle (see (3)), and the solutions hold for arbitrary SLE 
parameter k > 0. They appear in the form (9) 

C{z) = f{t^,m)G{^,m) , 

where the algebraic prefactor / is given by (10) and (11). It is independent of boundary conditions, and in fact a 
function of x and y via the coordinates ^ = ^{x, m) and i}} — ipiy, m), given by elliptic functions (see (8)). ^ and V' are 
the natural coordinates for this problem, as explained in [1]. The factor G is a single conformal block, given by an 
algebraic factor and an Appell hypergeometric function (see (12-16) and (18)). G depends on boundary conditions, 
and determining how the various Gs contribute to the cluster densities for the different boundary conditions of interest 
is one of our main results. (G can depend either on ^ and m, as written here, or alternatively on ip and m as mentioned 
in subsection II A, depending, respectively, on whether the right and left or top and bottom sides of the rectangle are 
'wired.') 

Subsections II B and II C next examine the question of expressing cluster densities (resp. crossing probabilities) for 
the boundary conditions of interest (sec figure 2) in terms of the Gs (resp. hypergeometric functions), and also explains 
how the various cases are implemented in specific models. Making use of results in [21], an explicit form for a horizontal 
crossing probability that generalizing Cardy's horizontal crossing probability for percolation to these models is given 
in (30) (see also the lines just below this equation). Expressions for the density Pi{z^R) (resp. Pr{z, R), Pir{z, R)) 
of clusters anchored to the left (resp. right, left and right) side(s) of a rectangle with either independently wired left 
and right sides in terms of the weights 11^, TlBg, 11^^, Tl^ and Ily (and the corresponding densities Pi{z,R) (resp. 
Pr{z, R), P^riz, R)) for mutually wired sides) are given in (31)-(33). (Explicit formulas for Hh and Hy are found in 
(54), (55), (60) and (61).) Subsection II C then employs Coulomb gas methods to find solutions for the weights 11^, 
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and Hb^ (the first row of figure 2). Tliese are given in the three equations (43)- (45), 



IIa = /(^,-0,™)Gvi(C, m) 
Hs,. = /(e,V,™)Gi(e,TO) 



and 



The expressions for the bottom row of figure 2 follow by symmetry. This completes our determination of the densities 
for the models and boundary conditions of interest here. 

Section III then discusses the factorization behavior implied by our solutions, which generalizes previous results for 
percolation to a variety of critical models. For percolation, the universal ratio (49) 



was considered [14], with Ph{R) the horizontal crossing probability given by Cardy. In that case, p is constant to 
within 3% inside the rectangle, being equal to 1 when z is on the left or right side, and rising to about 1.03 when z is far 
from the sides. The numerator, a six-point correlation function, therefore factorizes into the lower-order correlations 
in the denominator to very good approximation. In this section, we generalize p to an expression (50) defined in 
terms of correlation functions restricted to be single conformal blocks by specification of the propagating channels. 
This then leads to the expression for p (63) in terms of the weights 11^, Hs^,, Hs^, XIr and Ily or equivalently the 
conformal blocks G and the quantities Gi.i and Gi^a, valid for all the critical models mentioned. In all cases, by 
definition p = 1 when z is on the left or right side. It is also constant when z it is far from the sides, but how much 
that constant deviates from 1 depends on the model (more exactly, on k). Thus the factorization changes significantly 
with X in some models. However, in all cases p is independent of y, as in percolation, because the algebraic prefactor 
/, which also depends on y via ip, divides out of the ratio, leaving only conformal blocks G and other factors that are 
independent of y. This y-indcpendcnce is a consequence of the unusual symmetry of the conformal blocks G found in 

[!]• 

Next, section IV makes use of the locality property of percolation to derive the explicit result (78) giving the density 
of percolation clusters in a rectangle with free boundary conditions on all sides 
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with F\ the Appell function (17) and K'im) = K{1 — m), with K the complete elliptic integral. 

In section V we compare our predictions for the ratio p (see (63)), which governs factorization, for the Q — 1, 2 
and 3-state Potts models, including both FK and spin clusters, for various aspect ratios. The agreement is very good. 
The formula (78) for the density of crossing clusters in percolation is also verified. 

The Appendix presents forms for the densities and crossing weights in models where the hypergeometric functions 
simplify. 
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Appendix A: Special parameter values 



In this appendix we consider special values of the parameter k for which any of the hypergeometric series that occur 
in the density conformal blocks or crossing weights (see (12-16) and (54-55) respectively) may be written in terms 



22 



of simpler functions. In particular we look at k values for which the functions have non-positive integer parameters, 
as this may imply that the series truncates into hypergeometric series in fewer variables, potentially even into a 
polynomial. These five quantities are either proportional to the Appell function Fi , which appears in the forms 

. ^ f 4 4 „ 16 „ 8 

l--;-,2-— ;2-- 



s,t] (Al) 
s,t) , (A2) 



_ r(2-8/^)r(16^-l) ^ f^_^ i i ^ 
' r(12/«:)r(l - 4/«;) k'k'k'k 

in (13, 14, 15) and (12, 16) respectively; or as one of the functions Gi,i{m) and Gi^^^m), defined in (54, 55), which 
are given by a hypergeometric function 2^1- In what follows we consider values 2 < k < 8 for which either 2 — 8/k, 
2 - 12/k or 2 - 16/k e Z < 0. There are ten such values: n e {8, 6, 16/3, 4, 16/5, 3, 8/3, 16/7, 12/5, 2}, which 
correspond to n e {0, ±1, ±\/2, ±2} in the critical 0(n) loop model dilute phase and n G {0, 1, \/2} in the dense 
phase. 

The common pre-factor given in (10) is already an algebraic function in ^ and -0. Thus, for brevity we will not 
always include it, and restrict our attention to the functions (12-18) and (54, 55). 

We also note that in some cases hypergeometric scries with rational parameters of small denominator can be written 
in terms of other relatively simple functions. The ten k values that we discuss are all those with 2 < k < 8 that 
give rational hypergeometric parameters with denominator of four or less. Because the denominator value has a 
pronounced effect on the type of functions that appear upon simplifying the hypergeometric series, and because the 
associated values of n for each denominator of four or less all have the same magnitude, we group these values together. 



1. Integer hypergeometric parameters: n = 2 

a. Dilute n — 2: Gaussian Free Field 

Parameter value k = 4 corresponds to the critical Gaussian free field or to the Q = 4 state Potts model. For this 
value the various hypergeometric functions are 

/ N 2-s-t(2-t) , , 1 

ui{s,t) = 2(T~s) "2(s,t) = - 

r< I \ n ( \ 

GiAm) = Gi^m) ^ — . 

The conformal blocks are proportional to 

Gi(^, m) = -, Gn{t m) - (1 - m)(l -f mg^) 



^J^m{l - m)C(l - ^)(1 - m^) ' ^J 4m{l ~ m)^{l - ^)(1 - m^) 
Giii(C, m) = — = — = = Gw[^, m) = —r== == — == 



V4to(1 - ™)C(1 - 0(1 - mC) yjAm{l - m)i{l - - 

^ (, ^ m(l-0' „ . {l-m){l-me) / ~ 

y/4m{l~m)^{l-C){l~m^) V4to(1 - to)^(1 - ^(1 - V m sn {2xK \m) 

The expression (63) for p becomes 



I + m 

P^""' " V l + mcn2 {2xK'\m) ' ^ ^' 



where double argument identities for the Jacobi elliptic functions have been used. 



6. Dilute n = —2; Loop-Erased Random Walk 



An SLE with parameter k ~ 2 has been shown to correspond to the loop erased random walk [27]. However, 
the loop fugacity is less than zero for that model, so the present correlation function does not translate to a set of 
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densities. As k 2 the hypergeometric functions become 

, , l + s-3t (l-tf , , , , s + t-3 
Ms,t) = + j^^—j^il-2s + 2t-st) U2{s,t) = 

4 4 

Gi,i(m) =l-2m + m^ — Gi,3(m) = -m^ + — 

And we identify a few key functions: 

m(l - m)^^ [3 - ^ - m^] 



Gi(e,™) = 
Gv(^,m) = 
Gvi(?,to) = 



2[e(i-e)(i-mc)]'/2 

77i(l - C)M2 - m + ^ - 2to^] 



(1 - m) [1 + m - 3mC + 3771^^5 - m^^^ _ ^4^6j - m) [4 + 4m - 3msn2(2iir'x|m)] 

2m2[e(l-e)(l-"^0]^^^ ~ m2sn3(2/4''a;|m) 

2. Half integer hypergeometric parameters: n — 

For fugacity n = in either the dense (k = 8) or dilute (k = 8/3) phase p{x,R) = 1, as discussed. A pecuhar 
property of n = models is the absence of bulk loops. Therefore every point (except for those that belong to the 
SLE hulls) is adjacent to some part of the boundary, so that the bulk density operator is a generalization of the zero 
weight indicator operator used by Schramm in [28]. 

Schramm's paper considered a single SLE hull in a simply connected domain, and calculated the probability that 
the SLE passed to the left of a marked point in the bulk. This is equivalent to a correlation function with two SLE 
operators to generate the SLE hull, and a bulk indicator at the marked point that returns one if the hull passes 
to its left and otherwise. If it is possible (not possible) to draw a path from the point to the right boundary 
without crossing the hull the hull passes left (right) of the marked point. Thus, when n ~ and the only loop in the 
configuration is the hull itself the interpretation of the indicator and the density operator are identical, as evidenced 
by the weight 2/ii/2,o = 0, which is also the the weight of Schramm's operator. 

In the current correlation function, the bulk operator acts as an indicator operator in the presence of the two SLE 
hulls generated by the four 0i^2 operators. The indicator distinguishes between six cases: three correspond to the 
possible locations of a point relative to a horizontal crossing of hulls (configurations A, Bf,, and Bf in figure 2), and 
three correspond to the analogous case for vertical crossing (configurations A, Bi, and Br). 

The sum of the weights with horizontal crossings equals the total weight of horizontal crossings. A similar result 
holds for vertical crossings. These relationships can be written as 

cx Ha + Hb, + Hb, (A4) 
Hy cx n35- + nB^ +nB, . (A5) 

The weights from the bottom row of figure 2 follow from (43-45) upon letting (^, m. A, Br, Bg) -> (-0, 1 — A, Bt, B^). 
To emphasize the relations between the densities and crossing weights, we include the pre- factor (10) in this section. 



a. Dense n = 0: Peano Curve 



We first examine the dense phase with k = 8. This case is equivalent to a space-filling Peano curve. The relevant 
hypergeometric functions reduce to elliptic integrals on applying well-known identities: 



ui{x,y) = Gi,i(x) = Gi,3(a;) = 2F1 



1 1 



1 



'-K{x) 



U2{x,y) ^ -Fi 



11 13 

2' 2' 2' 2 



x,y 



F sin"^ ^/^ 



where F(-|m) is the incomplete elliptic integral of the first kind. This satisfies F (sin ^ sn(z|m)|m) = z, i.e. it acts as 
an inverse elliptic function. 
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For K = 8 the function fK=si£,, V': ^) = {K') ^[m{l — m)] -^^^ and we find: 
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TT 
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-F f sin ^ = — 



F sin 




yF (sin^i VV' 



F sin 



We simplify these expressions using elliptic function identities that are related to the symmetry operations in (22). 

The value k = 8 corresponds to a space filling Pcano curve, equivalent to the the boundary between the uniform 
spanning tree (UST) [27] anchored to the vertical sides and its dual, which is another UST anchored to the horizontal 
sides. As the curve is space filling, any infinitesimal neighborhood intersects both the UST and its dual, so the 
densities of spanning trees and dual spanning trees are equal and uniform throughout the rectangle. 

Only one of the UST and its dual can be a single component tree, since by connecting the left and right edges we 
prevent the connection of the top and bottom and vice versa. Thus when the UST spans the rectangle from left to 
right there are two components to the dual clusters, one attached to the top edge and the other to the bottom. Along 
with the uniform density of the trees, this means that Hh oc Ha = + IIb^, the same as (A4). A similar result 
holds when the left and right sides belong to different trees, which implies a vertical crossing by the dual tree. 



b. Dilute n = 0: Self Avoiding Walks 



The value of k = 8/3 corresponds to the problem of self avoiding walks. For k = 8/3 + e the hypergeometric 
expressions all diverge with r(2 — 8/k) ~ — 8/9e, as e — >■ 0. We adjust the normalization of wi,W2,G'i_i and G1.3, 
dividing by a factor of — r(2 — 8/k) to eliminate this divergence: 



Gia(m) := 
Gi,3(m) := 
ui(s,t) := 



Gi,i(m) _ 15 
-r(2-8/K) ~ 32 
Gi,3(to) 



7 3 
2' 2 



;3 



1 — m + 



E{m) 



2 — 3in + 
2^ 



K{m) 



-r(2-8/K) 

uiis,t) 
-T{2-8/k) 



Gi,i(m) 

{s-2st + tY 8t{s-t){l-t) d 



s2(l-s)2 



5s(l — s) ds 



Gi,i(s) , 



where E{m) is the complete elliptic integral of the second kind. Thus similar to k = 8, the gaussian hypergeometric 
functions arc replaced by complete elliptic integrals. 
The crossing weights are 



n^(m) = {K'f 
Hy = {K'f 



bj(l — m) K(m) 



1 — m + 



E{m) 



27r 

2 — 3m + 
'hi 



K{m) 



The density pre-factor is 



{K'f[m{l-m)f' 



ds''{2K'x\m) + ds^{2K'y\l - to) 
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which, along with the simphfied expressions for the various G(^, m) functions, gives us 

^ ^ 4ds^(2i^'a:|m)ny - (8/5)to(1 - m)5,nny 

^ ds^{2K'x\m) + ds^{2K'y\l ~ m) 

_ 4ds^(2iir'a:|m)n// - (8/5)m(l - m)dmTlH 

^ ~ ds2(2is:'x|m) + ds^(2/i:'y|l - m) 

^ ^ Ms\2K'y\l - m)ng + (8/5)m(l - m)a,„ng 

^' ds2(27^'a;|m) + ds^{2K'y\l - m) 

^_ _ 4ds^(2A"y|l - m)ny + (8/5)m(l - m)dm'nv 

^ ds2(2is:'x|m) + ds^(2X'y|l - m) ' 

after simphfying with double argument identities for the elliptic functions. With the densities written in this way, 
relations (A4) and (A5) follow immediately. 



3. Hypergeometric parameters with denominator three: |n| = 1 

For these particular parameters (k = 12/5, 3 or 6) the hypergeometric functions do not simplify beyond the 
truncation of Gi.i(to) that occurs because > 2 — 12/k e Z. 

The truncated functions, which are related to the identity channel of the crossing weights (see (59)) are, 

{1 K= 6 

(l-27n)(l + m)(l-m/2) k= 12/5 

The SLE parameter k = 6 corresponds to critical percolation, or equivalently the dense phase of the 0(1) loop gas. 
This is a purely probabilistic limit. Here the form of the identity channel crossing weight is such that the sum of all 
crossing probabilities is one. 

The SLE parameter k = 3 corresponds to spin clusters in the critical Ising model, or equivalently the dilute phase of 
the 0(1) loop gas. This is proportional to the Ising partition function in a rectangle with the horizontal and vertical 
edges fixed in opposite spins. 

The SLE parameter k = 12/5 corresponds to the dilute phase of the 0(— 1) loop gas. 



4. Hypergeometric parameters with denominator four: |ti| — y/2 

For K — 16/3 or 16/5, which correspond to = Q = 2, our expressions simplify greatly. 



a. Dense n = y/2: Ising FK clusters 

An SLE with parameter k = 16/3 corresponds to the boundary of the critical Q = 2 FK cluster model. For k = 16/3 
the hypergeometric functions become 



l + ^/l-m, l-y/l-m 
Gis(m) =\ Gi,3(m) =\ 



2(1 -s) y 2s(l-s) 
The functions that simplify become 



(1 - to)3/8 1 + 

V2mi/8 [^(1 - ^)(1 _ m^]'^^ v/TTT^ 

(1 - m)3/** 1 - 

V2mi/8 [^(1 _ ^t)(i _ ,71^]'/^ ^/^^^ 
Giii = [Gii + Giv] ■ 



Gii — 



Giv — 
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Combining these with 

Gvi = 



(1 - 7n)3/8 



gives the relatively simple expression 



+ 5r(i/4)2 ^U'i'i'i 



/l_F_yme _ _32e^ /m(l_K/m)7r A 3 3 9 
^^•^^"Vl-V^^ 5r(l/4)2\/ l-m^2 1 1^4' 4' 4' 4 



(A6) 



6. Dilute n = y/2 

An SLE with parameter k = 16/5 corresponds to the dilute phase of the O (\/2) loop gas. At 
hypergeometric functions become 



= 16/5 the 



ui{s,t) = - 



^ „3/2 Vl-Vi- 



Gi,i(m) = [2-3VT"^+2(l-m)] ^ ^ 

1 - VI 



3/2 



s (1 + V^) 
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16 mV8 [^(1 _ ^)(l _ mi)f''^ (1 - m)^/^ 
G„i = [Gii + Giv] /V2 . 



(A7) 

(A8) 
(A9) 



c. Dilute n = — -\/2 

The final case corresponds to k = 16/7, the dilute phase of the 0(— \/2) loop gas. At k = 16/7 the hypergeometric 
functions become 



Om(™)^-(=) 



5/2 4(1 - m)2 - 10(1 - to)3/2 + 11(1 - m) - 10(1 - mf/'^ + 4 
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In this case, the expressions being rather lengthy, we leave it to the interested reader to assemble the simplified 
functions and densities. 
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FIG. 9: Same as upper figure 6, except for Q = 3 FK cluster densities (k = 24/5). The simulations were performed on a rectangle 
of about 10^(x4 resp.) square lattice spacings for R = 1/2 (resp. 1, 2, 3). Here the plateau is at Cl'^.^^^/i^'* = 1.15470 (dashed 
line), from (68). 
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FIG. 10: Numerical results for p{x,R) for Q = 2 spin cluster densities (k = 3) against the CFT prediction (solid curve). The 
data points are simulation values of p{x,y,R) averaged over the y-coordinate with x,R fixed. Simulations were performed on 
a rectangle of about 10® square lattice spacings. Notice that both the theory and data plateau at C^'f.j^ 3 = 3.15123 (dashed 
line) given in (68). The R — 1/2 data are too small to show. This data is presented separately in figure (12). 




FIG. 11: Same as figure 10, except for Q — 3 spin cluster densities (k = 10/3). The simulations were performed on a rectangle 
of about 10®(x4 resp.) square lattice spacings for R = 1/2 (resp. 1, 2, 3). Here the plateau is at 2Clll-^ :^/{l + ^5) = 2.08229 
(dashed line) given in (68). See figure (12) for expanded R = 1/2 data. 
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FIG. 12: Results for p(x,R = 1/2) for Q = 2 (left) and Q = 3 (right) spin cluster densities. The CFT prediction is the solid 
curve. The Q = 2 (resp. 3) simulation was performed on a rectangle with area about equal to 10^(x4 resp.) square lattice 
spacings. 
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FIG. 13; Density of percolation crossing clusters at i? = 2. Solid curves are predictions from (78), data points from simulations. 
Data deviates more from the prediction near the center since Pa^ is almost flat there. See text in section V C for values. 



